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Abstract 

We investigate mixing of neutrinos in the i/MSM (neutrino Minimal Standard Model), 
which is the MSM extended by three right-handed neutrinos. Especially, we study ele- 
ments of the mixing matrix 0^/ between three left-handed neutrinos Va (a = e, (j., r) and 
two sterile neutrinos Ni [I = 2,3) which are responsible to the seesaw mechanism gen- 
erating the suppressed masses of active neutrinos as well as the generation of the baryon 
asymmetry of the universe (BAU). It is shown that Qei can be suppressed by many or- 
ders of magnitude compared with 0^/ and QtI, when the Chooz angle ^13 is large in the 
normal hierarchy of active neutrino masses. We then discuss the neutrinoless double beta 
decay in this framework by taking into account the contributions not only from active 
neutrinos but also from all the three sterile neutrinos. It is shown that N2 and give 
substantial, destructive contributions when their masses are smaller than a few 100 MeV, 
and as a results Qel receive no stringent constraint from the current bounds on such 
decay. Finally, we discuss the impacts of the obtained results on the direct searches of 
A^2,3 in meson decays for the case when iV2,3 are lighter than pion mass. We show that 
there exists the allowed region for with such small masses in the normal hierarchy 
case even if the current bound on the lifetimes of from the big bang nucleosynthesis 
is imposed. It is also pointed out that the direct search by using 7r+ — )• e+ + A^2,3 and 
— 7- + A'^2,3 might miss such A^2,3 since the branching ratios can be extremely small 
due to the cancellation in Qei, but the search by — )• /i"*" + A'^2,3 can cover the whole 
allowed region by improving the measurement of the branching ratio by a factor of 5. 



1 Introduction 



The extension by right-handed neutrinos is one of the most interesting physics beyond the 
Minimal Standard Model (MSM), since it gives a simple solution to the problem of the neu- 
trino masses confirmed by various oscillation experiments. Usually, right-handed neutrinos are 
introduced with superheavy Majorana masses and sizable Yukawa coupling constants in order 
to realize the seesaw mechanism [1], which accounts naturally for the smallness of neutrino 
masses. Furthermore, the decays of such right-handed neutrinos can be a source of the baryon 
asjmimetry of the universe (BAU) through the leptogenesis mechanism [21 13] . When the masses 
of right-handed neutrinos are hierarchical, the observed BAU requires the mass of the lightest 
one should be heavier than about 10^ GeV [4J. Although such singlet fermions provide simple 
and natural solution to the origins of neutrino masses and BAU at the same time, it is almost 
impossible to test them at experiments in the future. 

It should, however, be noted that right-handed neutrinos can bring about important physical 
phenomena, even when the scale of Majorana masses are so light to be produced in terrestrial 
experiments. One attractive possibility is the z/MSM |5l[6] in which three right-handed neutrinos 
are introduced with masses below the electroweak scale^J Interestingly, this simple model can 
explain the origins of neutrino masses, BAU and also dark matter of the universe at the same 
time. The Yukawa coupling constants are so small that the seesaw mechanism still works, and 
mass eigenstates of neutrinos are divided into two groups, active and sterile neutrinos. The 
flavour mixing of active neutrinos accounts for neutrino oscillations observed in experiments. 
On the other hand, three sterile neutrinos, A^^i, N21 and N^, solve cosmological problems in the 
MSM. 



One of the sterile neutrinos, A^^i, plays a role of dark matter [T3]cj When its mass is in 
the keV order, it can be produced by the so-called Dodelson-Widrow mechanism [T7], i.e., by 
the thermal scatterings through its mixing with the left-handed neutrinos. See Refs. |18] |27]. 
This dark-matter particle receives severe astrophysical constraints [28]. One important bound 
comes from the X-ray background [I9]-[2T] and the other comes from the cosmic structure 
at small scales like the Lyman a forest [29]-[32]. Even when these constraints are imposed, 
as shown in Ref. [27], the correct abundance of the dark matter can be obtained through the 
Dodelson-Widrow mechanism by invoking the resonant production [18] in the presence of the 
large lepton asymmetry. Notice that A'^i plays no significant role in the seesaw mechanism [3]. 
This is because Yukawa coupling constants of A'^i is so small that its contribution to the mass 
matrix of active neutrinos is negligible. Due to the very suppressed interaction the direct search 
of Ni at experiments is very difficult. However, it can be observed by the specific spectrum in 



The explanation of the LSND anomaly in this framework had been investigated in Ref. [7] . The non-minimal 
coupling of the Higgs field to gravity allows to realize the cosmic inflation p]. Further, various extensions of 
the model have been discussed [9|-[13j. 

"i^^It has been investigated various phenomenon in astrophysics by the sterile neutrino dark- matter, e.g., the 
explanation of the pulsar kick |15j. See, for example, a review of this issue |16j and references therein. 
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the X-ray background coming from the decay of A^i into active neutrino and photon. It has 
also been discussed the search in laboratory [33] . 

The rest two sterile neutrinos, N2 and N^, are responsible to generate not only the seesaw 
masses of active neutrinos but also BAU through the mechanism [31]. The flavour oscillation 
between N2 and in the early universe induces the separation of lepton asymmetry between 
left- and right-handed leptons and the asymmetry in the left-handed sector is partially con- 
verted into the baryon asymmetry by the sphaleron processes [3S]. The z/MSM realizes this 
baryogenesis scenario without conflict with the observational data of the neutrino oscillations 
when their masses are quasi-degenerate and in the range (9(O.1)-0(1O) GeV [6]. See also the 
recent analysis in Refs. [36l [37] . 

It is interesting to note that sterile neutrinos N2 and can be tested in various experiments 
as pointed out in Ref. [38|. This is crucially important to reveal the origins of the neutrino 
masses as well as the cosmic baryon asymmetry. For this purpose, we would like to study 
the mixing of sterile neutrinos with left-handed neutrinos {a = e, n, r) in this paper. The 
elements of such mixing matrix, B^/, are vital to discuss phenomenology of the i/MSM, since 
the strength of the interactions of sterile neutrinos is determined by them. Sterile neutrinos in 
the model possess the Yukawa interactions and also the weak gauge interactions via the above 
mixing after the electroweak symmetry breaking. Since the elements Qai are proportional to 
the Yukawa coupling constants Fai, both interactions are controlled by Qai- 

It should be noted that the mixing elements of and A^3 can take values varying by many 
orders of magnitude. This point had already been discussed by using the model with the 
lepton symmetry [32] • As we will show in Sec. [31 the mixing elements increase exponentially 
as Gq,/ oc exp(Imci;) for large Imu {u is a complex parameter in the neutrino Yukawa matrix) 
keeping the masses and mixing angles of active neutrinos unchanged. This enhancement leads 
to the various significant impacts on phenomenology of sterile neutrinos. For example, it gives 
the larger production/detection rates in the search experiments, the larger CP asymmetry in 
baryogenesis, and the shorter lifetime of A^2,3 making them cosmologically harmless. 

The purpose of this paper is, thus, to investigate the mixing elements Qai in detail, and to 
reveal how the elements of N2 and N3 depend on mass hierarchy, mixing angles and CP phases 
of active neutrinos in addition to the parameters of sterile neutrinos. Our analysis will show 
that there can be strong hierarchy among the mixing elements Gg/, G^/ and Qri depending on 
choice of the parameters. Especially, it will been pointed out that the strong suppression in the 
mixing elements of electron type Ge/ can happen for the normal hierarchy of active neutrino 
masses, which enlarges the allowed region of the model. Although we shall consider the z/MSM, 
the results of the mixing elements Qai in this paper can be applied to the general models of 
the seesaw mechanism with two right-handed neutrinos. 

Further, we would like to discuss the implications to two phenomena of sterile neutrinos. 
The first one is the neutrinoless double beta decay in which the mixing elements of active and 
sterile neutrinos play the crucial roles. This problem had already been discussed in Ref. [HI]. 
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We will extend the analysis especially when the masses of N2 and A^3 are smaller than a few 
100 MeV, and show that the rates of the decays in the i/MSM is smaller than those in the usual 
case where only active neutrinos give the contribution. Thus, the model is free from significant 
constraints discussed in Ref. [H]. The other one is the search of A'"2.3 produced in the decays 
of charged pions and/or kaons where the elements |6e/P or |0^/P determine the production 
rates. Especially, we will point out that N2 and which masses are smaller than pion mass 
are still allowed by the constraints from the direct searches as well as that from the big bang 
nucleosynthesis [121 SS] for the normal hierarchy case. This is different from a conclusion from 
Ref. [38]. The reason for it lies in the cancellation in the Be/ mentioned above. In addition, we 
shall present some implications to the future searches by using the decays of vr"*" and . 

This paper is organized as follows. In Sec. [2] we briefly review the framework of the present 
analysis, i.e., the i/MSM. In Sec. [3] we study the mixing of sterile neutrinos N2 and in the 
charged current interactions. Especially, we investigate how the mixing elements depend on the 
parameters of active neutrinos, i.e., the mass hierarchy, mixing angles and CP violating phases. 
We then apply the obtained results in phenomenology of N2 and N^,. In Sec. S] we estimate 
the contributions of sterile neutrinos to the neutrinoless double beta decay and address the 
importance of such contributions when the masses of N2 and A^s are smaller than about the 
order of 100 MeV. In Sec. [5] we discuss search of N2 and in the charged pion and kaon 
decays for the case when their masses are lighter than the pion mass. Finally, our results are 
summarized in Sec. [61 We add App. |X]to present the expressions for the mixing elements. 

2 The z/MSM 

First of all, we review the z/MSM, which is the MSM extended by three right-handed neutrinos 
'^R/ (-^ = 1) 2, 3) with Lagrangian 

C-vMSM = ^MSM Vri - (^Fal ^ UrI + ^ Vr^j Vrj + h.C^j , (1) 

where £msm is the MSM Lagrangian. $ and La = (eiayJ^La)^ (a = r) are Higgs and 
lepton weak-doublets, respectively. We denote Yukawa coupling constants of neutrinos by Fai. 
Here and hereafter we work in a basis in which the mass matrix of charged leptons is diagonal. 
In this model neutrinos receive the Majorana masses [Mm]ij = Mj6jj (which are taken to be 
real and positive without loss of generality) and the Dirac masses [M£)]ai = (($) is a 

vacuum expectation value of the Higgs field). The distinctive feature of the model is the region 
of the parameter space of Eq. i^, i.e., we restrict ourselves in the region 

|[M,j],/| <M, < Aew (2) 

Notice that the seesaw mechanism still works even if the Majorana masses are smaller than or 
comparable to the weak scale Aew = C'(IO^) GeV. This is simply because neutrino Yukawa 
coupling constants of interest are extremely small. (See the discussion below.) 
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The mass matrix of neutrinos M, which is a 6 x 6 symmetric matrix, is given by 

Md 
tl Mm 



^=1 Ml Ml )• 



We can diagonahze it by using the unitary matrix U as U'^ M U* = M^"^s_ ^i^e seesaw mecha- 
nism shows that If at the leading order takes the form 

.?V (4) 



-etf/ 1 ^ 

Here U is the 3x3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [H]; 

U'^MyU* = diag(mi, m2, ms) , (5) 

where M^, = —MdM^^MJ^ is the seesaw mass matrix. We call the eigenstates having masses 
rrii as active neutrinos Ui {i = 1,2,3). The rest three mass eigenstates, denoted by Nj, are 
almost corresponding to right-handed neutrinos Nj ~ urj having masses Mj. The neutrino 
mixing in the charge current is then induced through 

l^La = Uai + Qal Nj , (6) 

where the 3x3 mixing matrix B is found at the leading order as 

9 -IM^ (7) 
0.1 - -jj^ , (7) 

and hence |0a/| ^ 1 due to Eq. ([2]). We shall call Nf as sterile neutrinos since they possess very 
suppressed gauge interactions. It should be stressed that sterile neutrinos here are originated 
from right-handed neutrinos in the seesaw mechanism. 

In the z/MSM three right-handed neutrinos play important roles in cosmology. One of them, 
say A^^i, is a candidate for dark matter of the universe. This dark-matter particle receives severe 
astrophysical constraints as mentioned in Sec. [TJ Even then, the correct dark-matter abundance 
can be obtained through the mechanism [T7] with the resonant production [18] in the presence 
of the large lepton asymmetry. The recent study shows that the required mass is Mi = 4-50 
keV and the Yukawa coupling constants are typically \Fai\ = 5 x 10~^^-4 x 10"^'^. As a result, 
A^i gives no significant contribution to the seesaw mass matrix Mj, [5]. 

The other right-handed neutrinos, N2 and A^3, are then responsible to the masses and mixing 
of active neutrinos. Notice that in this case the mass of the lightest active neutrino becomes 
mi < C(10-^) eV. Further, N2 and N3 can explain the origin of BAU. The flavour oscillation 
between them in the early universe can generate BAU through the mechanism proposed in 
Ref. [31]. In the i/MSM the correct amount of BAU can be obtained when N2 and are 
quasi-degenerate in mass [51 [351 EZ] • 

The main purpose of this paper is to study the mixing elements Qai {I = 2,3) of sterile 
neutrinos N2 and N3 with flavour neutrinos. To do this, let us express their Yukawa coupling 
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constants by using mixing angles and masses of active neutrinos in oscillation experiments. As 
mentioned before, the successful dark matter scenario requires very small Yukawa couplings of 
A''! and its contribution to can be neglected. Thus, we set F^i = here for simplicity. (See, 
however, the discussion in Sec. HI) In this case the neutrino Yukawa matrix F for N2 and N^, 
which is a 3 X 2 matrix, can be expressed without loss of generality as [l5l |16] 

F = ^^UDl/^nDll\ (8) 

Here parameters of active neutrinos are their masses D^, = diag(mi, m2, ma) and the mixing 
matrix 

C12C13 S12C13 sise"*"^ 

U = I -C23S12 - S23Ci2Si3e*^ C23C12 - S23Si2Si3e*'' S23C13 | X diag(l , e*'' , 1) , (9) 

S23S12 - C23Cl2Sl3e*'' -S23C12 - C23Si2Si3e*'^ C23C13 

with Sij = sin 6ij and Cij = cos 6ij. Note that there is one Majorana phase rj in addition to 
Dirac phase 6 under the considering situation. Because we have set Fai = 0, masses of active 
neutrinos are 

^3 = matm > m,2 = mgoi > TTii =0 in the NH case , 

— V ^atm > ^3 — in the IH case , (10) 

The observational data of mixing angles are s^2 = 0.318;;'^o;Q4g, S23 = 0.50^0^4, and < 0.053, 
respectively, and masses are m^^y = Amg^ = (7.59^Q;5g) x 10^^ eV^ and m^^.^^^ = jAmg^^l = 
(2.40lo:D X 10"^ eV^ (at the 3a level) |17]. Hereafter, we shall adopt the central values unless 
otherwise stated. 

On the other hand, parameters of N2 and are their masses Dn = diag(M2, M3) and the 
3x2 matrix 



^7 = I cos tu — sin tu I in the NH case 
^ sin w ^ cos u 



cos u — sm u 

= \ .^sina; .^cosa; | in the IH case , (11) 


where ^ = ±1 and u is an arbitrary complex number. Notice that the change of the sign ^ can 
be compensated by w — —u together with the redefinition of as ^N^ ^ [IB]. 



3 Mixing matrix of sterile neutrinos 

The important parameters for phenomenology of sterile neutrinos A^2 and A^3 are their masses 
M2,3 and mixing matrix 0. Especially, the latter one is crucial to specify the strength of 
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interactions with other particles. Here we would like to discuss how they depends on the 
parameters of active neutrinos. 

Before discussing the z/MSM, let us consider a toy model with one pair of left- and right- 
handed neutrinos. In this case, the mixing of sterile neutrino is given by [cf. Eq. (jlj)] 



lei^ = = ^ = 4.9 X ^ , (12) 

' ' Mtv V y V2.4 X 10-3 eVV 

where Md and are the Dirac and Majorana masses, and we have used the seesaw formula 

for active neutrino mass M^, = \ — M1)/Mn\- Thus, the mixing is determined from the masses 

of active and sterile neutrinos. 

In the i/MSM, since the parameter space is larger, the mixing matrix of N2 and is more 

complicated and its elements can be much different from Eq. (IT^ . Especially, as pointed out 

in Ref. the larger mixing can be obtained in the model with U(l) symmetry. We shall 

reanalyze this point by using the parametrization of the Yukawa matrix F presented in Eq. ([8]). 

The key for this issue is the complex parameter u in the Q matrix. It can be seen that the 

Yukawa coupling constants as well as the elements of mixing matrix become exponentially large 

as Fai, Qai oc expdlmwl) for |Ima;| 1, as long as the seesaw approximation is valid. It should 

be noted that the tiny neutrino masses observed in the oscillation experiments can be obtained 

even in this case. To express this enhancement factor, we introduce a parameter X^^ by 

X^ = exp(Imw) . (13) 

Before going into details, let us here summarize the general properties of mixing elements 

(i) |0a/|^ can be divided into X^, X° and X~'^ terms. 

(ii) The X^ term in |0q,2|^^2 is exactly the same as that in IBasl^Ms for a = e, n, r. Similarly, 

9 I 1 2 I 1 2 

the X~ term in |0q,2| M2 is exactly the same as that in |0q3| M3. 

(iii) The X° term in |0q,2|^^2 is opposite to that in |Bq,3|^M3 for a = e, /i,r. 

(iv) The coefficient of the Xj^ term in |Ba/|^M/ is obtained from the X^ term by changing 
,^ — 7- — ,^ for a = e, /i, r and J = 2, 3. 

We have confirmed these properties by direct calculations. From now on, we will present the 

1 2 

expressions of \Qai\ for \lmuj\ ^ 1 and discuss how they depend on the neutrino parameters, 
namely mass hierarchy, mixing angles, and CP violating phases of active neutrinos. 

We first consider |9q,/|^ for X^^ ^ 1 {i.e., Imu ^ 1) in the NH case. The leading 0{X^) 
terms are found as 

cos^6'i3 [tan^6'i3 + 2A/r^^sin(5 + 77) sin6'i2 tan 6113 + sin^ ^12] ,(14) 
sin^ ^23 cos^ ^13 [1 + 0(v^)] , (15) 
cos^ ^^23 cos^ [1 + 0(v^)] , (16) 



|0e/|'| 




= xl 


^atm 
4M7 






= xl 


4M7 






= xl 


^atm 

4M/ 
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where = '^soi/'^atm — 0.18. (The complete expressions for the 0{X^) terms are collected 
in App. |Al) Notice that these expressions hold for both 1 = 2 and 3 thanks to the general 
property (ii). It is seen that all the elements are proportional to X^^matm/^/, and hence they 
can be much larger than Eq. f|T2|) for ^ 1. 

Since the experiments show that is small and 623 is close to 7r/4, 16/^/1^ and |0r/|^ can 
be determined as 



|e,,|' ^|0.,|^^x^!^ = 6.1 X 10-1^ (^^] . 

' ^ ' ' ' 8M1 \ Mj J 



(17) 



On the other hand, the element Qei behaves quite differently. Indeed, it is interesting to note 
that the X^, terms in |0e2|^ and [Oesl^ vanish at the same time, when 



esin(5 + r/) = -l, (18) 
and the mixing angle 6^13 takes its critical value 

tan 9^1 = ^/r^sm6l2 ■ (19) 
The experimental data of 9i2, "n^atm and msoi with 3a errors gives the critical value of 613 as 



sin = 0.041-0.070, which can be below the 3a upper bound sin ^13 < 0.053E3 Moreover, 
we find in this case that the 0{X^) term also vanishes and only the (9(Xj^) term is left as 

\ ell ^ Mi \ Mi ) ' ^ ^ 

which becomes much smaller than other elements in Eq. (fT71) for ^ 1. We should mentioned 
that the above cancellation in |0e/|^ can be realized for any choice of masses M2 and M3. As 
we will show in Fig. [H the strong suppression in \Qeif is still possible when 6^13 is close to 611- 
On the other hand, when 6*13 = 0, |9e/| receives no strong suppression described above, but it 
satisfies the relation 

2r^sin2^i2 ~ 0.11. (21) 



This relation had already been obtained in Ref. [26j. 

In Fig. [T]we show the mixing elements |0q,2|^ in terms of X^ when M2 = 120 MeV. First, we 
observe that |0^2|^ and |6t-2|^ scales as X^ for X^ ^ 1, and they can take much larger values 
than the naive result in Eq. f[T^ . This behavior does not change much as long as 613 lies in the 

9 I 1 2 I 1 2 

experimentally allowed region. Second, when sin 613 = 0, |9e2| behaves similar to |6^2| and 
I Gt-2 I ^ for X^ ^ 1 , but is smaller by one order of magnitude as shown in Eq. fl^Tl) . Finally, it is 
clearly seen that |6e2| can be suppressed by many orders of magnitude for X^^ ^ 1 when ^^13 
becomes close to its critical value with a suitable parameter choice. Moreover, if 6*13 = 6^^^, we 



#3when we use the data at 2a level, sin^ 9% = 0.046-0.065, which exceeds the 2a bound sin^ 6*13 < 0.039. 
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Figure 1: Mixing elements |0o2p in the normal hierarchy in terms of X^^. We take sin^ ^^13 = 
(left), 0.05 (center) and siv? (right), respectively. The red solid, green dashed and blue 
dotted lines correspond to |0e2|^5 I^msI^ ^^"^ |0r2|^5 respectively. Here we take M2 = 120 MeV, 
Rew = 7r/4, 5 = n/2, rj = n, and ^ = +1. 



can see that \Qe2f is proportional to X^'^ and it can be extremely suppressed for X^ ^ 1. The 
cancellation of \Qei\ in the NH case is one of the most important observation in this analysis 
and we will discuss its impacts on the experimental signatures of N2 and later. 

On the other hand, when <^ 1 {i.e., Imoj <^ "1), the leading order contribution is 
proportional to and their expressions are given by Eqs. fll4l) . (fT5i) and (|T6l) by replacing X^ 
by X'"^ and ^ by — ^ due to the general property (iv). Therefore, the above arguments with 
the opposite sign of ^ exactly hold and we will discard this case in the followings. 

Next, we turn to discuss the IH case when ^ 1. In this case, the leading term of the 
mixing element of electron type is found as 



19 



el\ 



^>2 

4Mi 



COS^ 012 COS^ 6^13 



tan^ 012 — 2^ sin 77 




mi mi 

tant^i2 H 

m2 



(22) 



and the expressions for |0^/|^ and |0r/|^ are so long and they are collected in App. |X1 It is 
interesting to note that the X^ term as well as the X° term in |0e/|^ vanishes when 



^ sin?7 = +1 and tan6'i2 = tan^ 



12 




r 



-1/4 



(23) 



and then 6*12 is close to the maximal angle 7r/4. Unfortunately, it is far beyond the current 

2 

data of 3(7 range, and the cancellation in \Qei\ cannot be realized in the IH case, which is 



This point gives significant effects on the discussions given in 



different from the NH casel*^ 
the following sections. 

However, sin 17" plays a crucial role to determine the mixing element |0e/|. To see this 
point, let us take 6'23 = vr/4 and 6'i3 = for simplicity. In this case, the mixing elements are 



When ^sinry = —1, ia.n9i2 — yjvaijvnx, and sin^ia 
the same time. However, the required value of Qyi is not allowed by the current data. 



0, the X"^ terms in both |0^/|^ and |Br/|^ vanish at 
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Figure 2: Mixing elements |0a2p in the inverted hierarchy in terms of X^. We take sin 77 = 
7r/2 (left) and 37r/2 (right), respectively. The red solid, green dashed and blue dotted hues 
correspond to |0e2|^5 [©msI^ ^^"^ |0r2|^5 respectively. Here we take M2 = 120 MeV, Reo; = 7r/4, 
5 = 7r/2, = 0, and ^ = +1. 



given by 



I Q^j I ' I ^ ( 1 - e sin sin 2^^i2 ) , (24) 

\QMx. ^ \Qri\\, ^ + e sin r/ sin 2^,2) , (25) 

" oik// 



which gives the relation 



\Qeif _ \<S>eif _ ^ 1 -^sinr/sin2^i2 _ / 0.071 for^sinr7 = +l 



,0/./!^ |0r7|^ l + ^sini] sin 2^12 1 56 for^sinr/ = -l 

Therefore, the mixing element of electron type can be smaller or larger than others depending 
on the choice of sin 77". This property is represented in Fig. O It is seen that |0/i/|^ and 

I 1 2 I 1 2 

10,-/! are almost the same, but |0e/| can be different from others. 

Before closing this section, we would like to stress again that the above results of the 
hierarchy between the mixing elements 0^/ are independent on the masses of sterile neutrinos. 
Therefore, they can be applied to the general seesaw models with two right-handed neutrinos. 



4 Neutrinoless Double Beta Decay 

The neutrinoless double beta (Oz/2/3) decay is one important phenomenon in which the mixing 
of active and sterile neutrinos, Uai and 0^/, plays a crucial role. The 0^2(3 decay in the i/MSM 
had already been investigated in Ref. [10], in which the contributions from active neutrinos 
and dark-matter sterile neutrino A^^i are estimated and those from A^2,3 are neglected since their 
masses are assumed to be so heavy that they decouple from the considering decay processes. 
Further, it had been discussed in Ref. ^1] that the Oz/2/3 decay gives the stringent constraint on 
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the mixing element Qei especially when their masses are around 100 MeV. Following to these 
analyses, we would like to reconsider this issue in this section. In particular, we shall take into 
account the contributions from N2^3 by considering a wider range of their masses, and also the 
mixing elements Oq2,3 discussed in the previous section. We will then show that the constraints 
from the Oi/2/3 decay are negligible in the i/MSM. 

The rate of the Oz/2/3 decay is characterized by the effective neutrino mass iries (see, e.g., 
Ref. gS]). In the z/MSM it is given by 

m,s = m:^+ J2 MiQljfp{Mi), (27) 

7=1,2,3 

where the first term m^^ denotes the contribution from active neutrinos 

<fr = E ■ (28) 

1=1,2,3 

The second term in Eq. f l7r|) denotes the contribution from sterile neutrinos in which we have 
introduced the function fp to represent the suppression of the nuclear matrix element from 
neutrinos with masses heavier than about 100 MeV [IHl ED]- In this analysis, for simplicity, we 
shall assume fp{Mi) = 1 for Mj < kp and fp{Mi) = {kp/MiY for Mi > where the typical 
energy scale in the matrix element is taken as = 100 MeV. The more precise treatment of 
the function fp at M/ ~ does not alter our final conclusions. 

First, we consider the case when all the sterile neutrinos in the i/MSM are lighter than A^. 
As for the dark-matter sterile neutrino A^i, its mass is indeed smaller than A^. On the other 
hand, the masses of A''2 and A^3 can be smaller or larger than A^. Here we consider the former 
case and the latter case will be discussed separately below. In this case the effective neutrino 
mass exactly vanishes. This is because 

meff= Yl MiQlj=[UM^'''^U^l, = [M,],, = 0. (29) 

1=1,2,3 /=1,2,3 

This cancellation has been recently observed in Ref. [50] by using the general seesaw model. 
Thus, the i/MSM with Mi^2,3 < kjs predicts zero event in the Oi/2/3 decay, which means that 
the current experimental limits on rries give no constraint on the model. 

On the other hand, A^2,3 can be heavier than kjs in the i/MSM. In this case, the prediction of 
meff is modified. We consider the case when A^^i is resonantly produced via thermal scatterings in 
the presence of the lepton asymmetries as mentioned in Sec. 2. In this already pointed 

out in Ref. [10], the contribution from dark-matter sterile neutrino A^i, denoted by m^g, is 
very small. The correct abundance of dark matter can be obtained when Mi = 4-50 keV and 
\F^i\ ~ 5 X 10-^^-4 X 10-13 [27]. since Mi < A^, we find 

= Ml \eli\ = ^-^^P^ = O{10-''-10-') eV, (30) 
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which is so small to be negligible in most cases, but we will point out its importance. 

Let us turn to consider the contribution from N2 and N^^. These sterile neutrinos are quasi- 
degenerate in order for the successful baryogenesis [Ml El [36l [37] . We then write the masses as 
M3 = Mat + AM/2 and M2 = Mn - AM/2 with AM < M^v, and divide m^^ from A^2,3 into 
two parts 

1=2,3 



J2 M.QlMM,) = mfj'^ + 6m^r , (31) 



where 



- MMn) J2 ^^^el , (32) 

1=2,3 



Sm^r = E M^^) - //^(^^)] ^^©e/ • (33) 
1=2,3 

Notice that the second part Sm^g '* vanishes when AM = 0, and more interestingly, that the 
first term can be written as 

m^r = -MM^)<E- (34) 
Therefore, the effective neutrino mass in the z/MSM when M2,3 > is written as 

meff = [1 - MM^)] m^ff + m^i + 6m^i'' . (35) 
It is then found that when Mjy ^ (and hence //3(M2,3) ^ 1) 

mcff ~ m^^f^ + m^^ , (36) 

which shows the sizable contributions come only from active neutrinos and dark-matter sterile 
neutrino |40]. On the other hand, A^2 and N3 give a significant effect on rries especially when 
they are quasi-degenerate and Mjy ~ A^, namely, they induce the destructive contribution 
to rries given in Eq. fl27j) . Thus, |meff| in the i/MSM can be much smaller than {m'^gl when 
Mjy ~ A/3. This point is illustrated in Fig. [3l where we show the allowed region of |mefr| in 
terms of Mjy by varying Recu, 6 and rj in the range [0,27r], Imo; in the range [0,7], and the 
parameters of active neutrinos within the experimental 3a range. 

Next, we turn to see how rrics depends on masses, mixing angles and CP violating phases 
of active neutrinos. In the NH case, the contribution from active neutrinos is {e.g., see [51] ) 

<ff = e-''^ (sin2 m3 + 6^'^^+"^^ cos^ ^13 sin^ 612 ^2) . (37) 

Notice that mi = in the limit of F^i — )■ 0. We also find that 

sm^r = ^f,{M^) [M3 QI3 - M2 ey 



.AMAl 



2M3 



^2 g-2iRe^ (^/^e sin ^13 - ze*(^+'')v^cos^i3 sin ^12) 



+ X~^e+^'^''^ (y^esin^i3 + ie^(^+'')V^cos^i3sin^i2)^ , (38) 
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Figure 3: Allowed regions for |meff| as a function of the mass of sterile neutrino M^, for NH 
(left) and IH (right) cases, respectively. The light gray region with solid lines, gray region with 
dot-dashed lines, and dark gray region with dotted lines correspond to regions allowed by the 
observational data for active neutrinos at the 3a, 2a and la level, respectively. Here we take 
AM/Mn = 5 X 10-*^. 
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where we have neglected the higher order terms of AM. 

It should be noted that m^g and also the term in 5m^'^ vanish, if cos 2(5 + t]) 
and tan^ 6'i3 = rmsin^6'i2. Interestingly, the required conditions are the same as f|T8|) and f|T9|) . 

eff 



i.e., those for the cancellation in |0e/p. In this case, Sm'^g^ becomes suppressed for a large 
together with a small mass difference of A'"2,3. As a result, the dark- matter sterile neutrino gives 



the dominant contribution to the effective mass in the Oi/2/3 decay |meff| ~ l^gfr I? which is very 
small as shown in Eq. (I5U]) . This is the reason why the 3a lower bound on |meff| is beyond the 
region of the plot. 

When the conditions (fT8|) and (fT9l) are not satisfied, we find for large 



I r -^2,3 I 



1.5 X 10"^^ eV 



(AM/Mn\ [I GeV\ 



N 



X' 



#5 



(39) 
This shows 



5x10-8^ 

by using the central values of parameters of active neutrinos and 5 + 7] = 7r/2, 
that we may neglect 5m^g^ as long as AM is sufficiently small. Note that X^ cannot be so 
large due to the experimental upper bounds on Qai, which will be investigated in the next 
section. In this case, the effective neutrino mass becomes 



m,s ^ [1 - fi^iM^)] + m^^ . 



(40) 



We then understand the 3a upper bound on |mefr| in Fig. [3] as [1 — /^(^at)] '^egljyj^x' "W'hich 
becomes suppressed when Mat is close to A/j. 
In the IH case, we obtain 



m 



Here the maximal value of \6m, 
prefactor becomes 1.8 x 10"^'^. 



- cos^ 6*13 (cos^ 6*12 mi + e^''' sin^ 6*12 7112) , (41) 
is shown. When S + rj = 3tt/2, it takes the minimal value and the 



N2.: 

off 



12 



and 



2 

AMAl 



2M% 



cos 6^13 



Xl e-2^^e^ {./^li cos + i e^^y^sin 612) ' 



+ Xj2e+2*Rea. (y^^ cos ^12 -ze^^yW sin ^12)^1 . (42) 



Then, as in the NH case, m^^ as well as the term in vanish, when ^ sin 17 = +1 and 

tan^ ^12 = mi/m2. The required conditions are the same as those of the cancellation in |0e/|, 
and they are not satisfied the present experimental data on the masses and mixing angles of 
active neutrinos. This is the reason why |meff | in the IH case receives the stringent lower bound 
as shown in Fig. [3l and m^^ becomes negligible. Further, compared with other terms, Sm^g'^ 
can be neglected for a small mass difference of A'"2,3. This is because 



when we take rj = 37r/2. |^ Therefore, we obtain in the IH case 

meff^ [l-/^(Mjv)]m:g. (44) 

In summary, we have shown that the effective neutrino mass in the 0i'2f3 decay is given by 
the contribution from active neutrinos m'^g (together with a very small m^i)forM2,3 >A/3. On 
the other hand, when M2,3 ~ A^, sterile neutrinos N2 and give the destructive contribution 
and then rries becomes much smaller than that from active neutrinos. Especially, meg vanishes 
when all the sterile neutrinos are lighter than A^. Therefore, the i/MSM receives no stringent 
constraint on the mixing elements Qai from the decay pointed out in Ref. [ilj . 



5 Search for light sterile neutrinos 

In this section we shall discuss the experimental search of sterile neutrinos A'"2,3. In particular, 
we restrict ourselves here to the case when the masses are smaller than the pion mass. This is 
simply because we would like to discuss the production of A^2,3 in the decays of 7r+ and K^. 
In addition, although such light sterile neutrinos have a long lifetime and may decay after the 
big bang nucleosynthesis (BBN) starts as we will show below, it can be avoided the stringent 
constraint from the BBN [52j on the hadronic decays of A'"2,3. In this section, we shall explore 
in detail the impacts of the results in Sec. [3] on the direct searches of A'"2,3 in the vr"*" and 
decays. 

"^^Here the maximal value of |i5m^'^| is shown. When 77 = 7r/2, it takes the minimal value and the prefactor 
becomes 8.0 x 10~^^. 
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Figure 4: Upper bounds on |Oe/P (left), |0e/0/i/| (center) and |6/x/P (right) from direct 
searches of sterile neutrino. In the left panel the red sold line is from [5l], the blue dotted line 
is from [55], and the cyan dot-dashed line is from [59]. In the center panel the red solid line 
is from [5^. In the right panel the red solid line is from [56], and the blue dashed line is from 

There have been so far various experiments of the direct search for sterile neutrino which 
give the upper bounds on the mixing elements Qai- In the considering mass region, the signal 
of sterile neutrino can be investigated by the peak search [53j in the energy spectrum of charged 
leptons from meson decays [5l|-[57|, and also by finding charged leptons from the decays of 
sterile neutrino inside the detector [58|-[60]. See, for example, Refs. [HI [381 [62]. The upper 
bounds on the mixing elements are summarized in Fig. [H Notice that the upper bound on 
|0r/| [60] is weaker than that on |0^/|, and then it plays no significant role for our conclusion. 
It can be seen that the bound on \Qei\ is severer than |0^/| by orders of magnitude. 

The mixing elements of A^2,3 scale as |0q,/| oc X^^ for X^^ ^ 1, as described in Sec. [3l Thus, 
the experimental upper bounds on |0q,/| can be translated into the upper bound on X^. Such a 
bound is the basis of finding the shortest lifetimes of N2 and N3 and also the largest production 
rates in the decays of vr"*" and K~^. Let us first discuss the case when Mn = 120 MeV. In 
this case, the experimental bounds are lOe/luB = 6.0 x 10~®, |6e7-6^/|uB = 1-8 x 10~^, and 

\e^i\l^ = 2.8 X 10-6. 

In the NH case, we show in Fig. \5\ the upper bounds on in terms of sin^ 6*13 since the 
element Gg/ is crucially dependent on this mixing angle. When there is no cancellation in Be/ 
(see the left panel of Fig. [5|), the experimental bound on \Qeif determines the upper bound 
on Xi^ in most cases. Note that, when ^13 is close to zero, |6e/P is smaller than |6^/P by one 
order of magnitude as shown in Eq. fl2Tl) and then |9e/6^/|uB puts the upper bound on X^^. 

On the other hand, as shown in the right panel of Fig. [5l the situation is drastically changed 
when the cancellation in Qei can happen as shown in Eqs. f|T8l) and f|T9l) . In this case, the 

I 1 2 I 1 2 

experimental bounds on \Qei\ and also |Oe/0^/| play no significant roles, but that on 
determines the upper bound on X^. Since the present limit |0/^/|ub weaker than others, 
the allowed region of becomes wider. Furthermore, we find that the CP violating phases 
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Figure 5: Allowed region in the sin^ plane for the NH case by taking S = and tj = ^ 

(left) and 6 = and ght). The red solid, the magenta dashed, and the blue dot-dashed 

lines are the upper bounds on from the experimental limits on |0e7p, |©e70/i/|, and lO/^/P 
respectively. The green line is the Sa limit on sin^^is. The black dotted lines show lifetimes of 
sterile neutrino. We take M3 = 120MeV, AMVM| = 10-^ Recu = f , and ^ = +1. 




Figure 6: Allowed region in the sin^ plane for the IH case by taking 5 = f and r] = ^ 

(left) and 5 = f and = f (right). The red solid, the magenta dashed, and the blue dot- 
dashed lines are the upper bounds on from the experimental limits on |0e/p, \QeIQ^lI\, and 
|0^/p respectively. The green lines shows the 3a range of sin^ O^. The black dotted lines show 
hfetimes of sterile neutrino. We take M3 = 120MeV, AM'^/Ml = IQ-^, 6*13 = 0, Reu = \ and 



change importantly the prediction of even when the condition f|T8l) is satisfied. This is 
because, as found from Eq. (!52|) in App. |Al the negative ^ sin 77 decreases |6/x/|- For instance, 
if we compare two sets of the CP phases, (5 = tt, 77 = 7r/2) and (5 = 0, 77 = 37r/2), satisfying 



the condition (ITSl) with ^ = the 
can be about 350 as shown in Fig. [Sl^^ 



atter case gives the suppressed for ^ 1 and X^ 



'^^''Here we have used the central values of parameters of active neutrinos. When we vary them in the Scr 
range, can be large as about 550 as shown in Fig. [7] 
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Figure 7: Allowed region in the Mj^f-X^^ plane for the NH (left) and IH (right) cases. The gray 
region is excluded by experiments for direct search of sterile neutrinos. The dotted lines show 
the lower bounds on for given lifetime of sterile neutrino. We use the observational data 
for active neutrinos at the 3a level. 

In the IH case, the mixing angle 612 is an important parameter to determine Qai and then 
we represent in Fig. [6] the upper bounds on X^ in terms of sin^6'i2. Since the cancellation in 

I 1 2 

Ge/ cannot be realized in the IH case, the experimental bounds on \Qei\ and |0e/0/x7| place 
the upper bound on X^. As explained in Eq. the choice of sin 77" is significant to 

determine the ratio |9e/|/|0;x7-|. When ^sinr] = —1 {e.g., ^ = +1 and r] = 37r/2), |6e7-| becomes 
larger than In this case, the upper bound on X^ becomes smaller due to the stringent 

experimental bound on |0e/| (see the left panel of Fig. [6]). Inversely, when ^ sin 17 = +1 {e.g., 
^ = +1 and T] = tt/2), \Qei \ becomes smaller than Then, the upper bound on becomes 

relaxed and X^^ can be large as about 100 (see the right panel of Fig. [6]). 

We then find numerically the upper bound on for a given M2,3 by varying Rew, S and rj in 
the range [0, 2tt] and also by varying the parameters of active neutrinos within the experimental 
3(7 range. The obtained results are shown in Fig. [71 In the NH case, the upper bound becomes 
weaker for lighter M^v. Especially, no sensible bound is available for M^r ^ 55 MeV. This is 
because the cancellation in 6e/ is possible and also because there is no stringent experimental 
bound on |6^/|^ for such mass region. It is quite interesting to note that X^ at C(IO^-IO^) is 
experimentally allowed for ~ 100 MeV. In the IH case, on the other hand, we can see that 
the bound on X^ is severer than that in the NH case. This is because there is no cancellation 
in 9e/ and X^ always receives the stringent experimental bounds on \Qeif and/or |0e/0^/| as 
explained above. 

The upper bound of X^ allows us to estimate the possible range of the lifetimes of N2 and 
N3. The lifetime is a key parameter from the following two reasons. One reason is concerned 
with the search for N2 and N^. In the test of sterile neutrino by finding charged leptons from its 
decay inside the detector, the lifetime determines the decay length and also the detection rate. 
The other comes from cosmology. The light sterile neutrinos under consideration can be long- 
lived particles and their decays would spoil the success of the standard cosmology, especially 



16 




20 30 40 50 60 70 80 90 100 110 120 130 20 30 40 50 60 70 80 90 100 110 120 130 

IVlN[MeV] MM[MeV] 



Figure 8: Lower bounds on lifetime of sterile neutrino in terms of Mat is shown by the red solid 
lines for the NH (left) and IH (right) cases. We use the upper bound on the mixing elements for 
sterile neutrinos, |6e/luB) l©M-fluB; |0e/0Ai/|uB and the observational data for active neutrinos 
at the 3a level. The BBN upper bounds on lifetime in Eq. fl46p are also shown by the blue 
dotted hues. 



the prediction of light elements of BBN. The lifetime is then restricted by various cosmological 
constraints. We shall discuss this issue later. 

In the considering mass region, N2^3 decay in to ViVjVj or {ia,i3 = e, /i) and their 

CP conjugate states. The lifetime, which is determined from the mixing elements B^/ and 
the mass of sterile neutrino, can be estimated by using the partial decay rates presented in 
Refs. [nH SSI ES] ■ Before dealing with the actual model, we consider the toy model described 
in Eq. f ll2p . In this case, since the dominant channel is the decay into three active neutrinos, 
the lifetime is estimated as 

~ ^92vr3 _ / rn^\ / lOOMeVA ' 

Therefore, sterile neutrino decays after BBN starts, which might lead to cosmological difficulty. 
In the i/MSM, however, this is not always the case and the much shorter lifetime can be achieved, 
because the lifetime scales as t^^ ,^ oc for ^ 1. Thus, a sufficiently large makes A^2,3 
cosmologically harmless. 

We find, as general properties, that the lifetimes of N2 and N3 are almost the same, and 
also that the lifetime in the IH case is slightly shorter than the NH case for a given X^. We 
then show the contour lines of lifetime of A^2,3 in the parameter regions discussed previously 
(see Figs. [5l [6] and [7]). It is seen that the lifetime can vary orders of magnitude depending on 
the choice of X^. Since we have obtained the upper bound on X^^ to be consistent with the 
search experiments, we can estimate the lowest value of the lifetime of A'"2,3. The results are 
shown in Fig. [H] for the NH and IH cases, respectively. In the NH case, the direct searches give 
no significant constraint on X^ for < 55 MeV. We then impose \Qal\ < 10-^ to ensure the 
validity of the seesaw mechanism. It is found that tatj 3 of C(10~^-10~^) sec is available in the 
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wide mass range. When the cancellation in Q^i occurs, the upper bound on gets larger. As 
a result, and |0t/| can take the larger values, which leads to the shorter lifetime of A^2,3 

even if \Qei\ is very small. In the IH case, on the other hand, the absence of such a cancellation 
gives the stringent lower bound of tatj 3 as 10-10"^ sec for Mn = 20-130 MeV. 

It had been pointed out in Ref. [38] that the bound from the BBN places the upper bound 
on the lifetime of A'"2,3, which results in the admitted window for the parameter space of the 
model. The decays of sterile neutrinos around the BBN era would alter the abundances of 
light elements in a desperate way. There are two significant effects by such decays. One effect 
is that the additional energy carried by sterile neutrinos increases the expansion rate of the 
universe. The other is the modification of the proton-neutron conversion rate. This is because 
active neutrinos produced by the decays of iV2,3 cannot be fully thermalized if the decays occur 
too late and then the distortion of the distribution functions affects the conversion rate. These 
two effects would alter the abundance of light elements (^He especially) too much, which leads 
to the upper bound on r^Vag. Notice that the stringent constraint from the BBN [S2j on the 
hadronic decays can be avoided in the considering situation. 

Such a bound had already been discussed in Refs. |l2l |33] , in which the authors consider the 
model with one flavour of active and sterile neutrinos. At present there is no analysis dealing 
with the realistic situation of the considering model. We shall postpone to study the issue 
in future publication [63], and just apply the bound on the lifetime given in Ref. [13]. The 
maximum allowed lifetime tbbn is given by 



where ti = 128.7, ^2 = 0.04179, and /3 = —1.828 [13]. (In this analysis we have taken a 
conservative bound for the mixing with z/^^^, and see the detail in Ref. j43j.) When = 120 
MeV, we find tbbn = 6.2 x 10~^ sec. 

In Fig. [8] the BBN bound on the lifetime tbbn is also shown in terms of Mn- We can find 
that, although such a short lifetime cannot be obtained in the IH case, there exits the allowed 
region in the NH case, which is different from the conclusion in Ref. [3S]. This is because of the 
cancellation in Q^i we have obtained in Sec. [31 which enlarges the allowed parameter space of 
the model by relaxing the stringent constraints on |0e/| from direct searches. In this allowed 
region, we can obtain remarkable predictions of the model parameters. This point is illustrated 
in Fig. [9] when = 120 MeV. First, should be sufficiently large for tat^ 3 < tbbn and 
Xi^ > 250 in this case. Second, in order to sufficiently suppress |0e/|, the parameters should 
be very close to those in the conditions ([T8|) and f|T9|) . Namely, the mixing angle di^ should be 
large (close to 6*^3), and the CP violating phases should be ^ sm{6 + ri) ~ 1. Finally, the negative 
sign of ^ sin?7 is required to decrease Therefore, the range of the Chooz angle 6*13 and 

the CP phases 6 and 1] is highly restricted. This leads to an important effect on baryogenesis 
of the model, which will be discussed elsewhere [63] . 





(46) 
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Figure 9: Allowed regions in the sin^ Oi^-X^^ plane (left) and 6-1] plane (right) when tatj 3 < tbbn 
for the NH case. Here we take M^r = 120 MeV. In the right-panel red points and green points 
correspond to the case ^ = +1 and ^ = — 1, respectively. 



In the following, however, we shall adopt the upper bound on the lifetime as tn2^^ < 1 sec, 
in addition to TTVag < tbbn; as the most conservative case, and consider the both NH and IH 
cases. 

Now, we are at the position to discuss the experimental search of A''2,3 by the peak searches 
in the decays vr"'' — > e"*" + A'"2,3, — >■ e"*" + iV2,3, and — > /x'*" + N2^3. The branching ratio for 
the production process is given by [53] 



BR(P+ + Nj) = BR(P+ + i^e) x K{mp, m^, Mj) 



where P = tt,K and 



e, /i. The function K is 



K{mp,mi,Mi) = \Q 



m 



2\2 



ei\ 



f3u 



22 4 
nipmi — m| 



where 



/3 



Ni 



rrii 



m. 



/3. 



1 - 



mo 



m. 



(47) 



(48) 



(49) 



Here we have neglected the masses of active neutrinos. It is then seen that the branching ratios 
are proportional to the mixing element squared, and hence X^. 

By taking = 120 MeV, the allowed regions of the branching ratios are shown in Figs. [10] 
and [11] for the NH and IH cases, respectively. Here we have imposed the lifetime bound 
in addition to the experimental bounds on the mixing elements, and the parameters of the 
Yukawa matrix have been varied as we did before. 

In the NH case, it is found that the branching ratios of 7r+ — )■ e+ + and — > e+ + can 
be very small even if one imposes the lifetime bound. This is due to the cancellation in |0e/| 
pointed out in this paper. As explained in Eq. (!20|) . |Oe/P is proportional to Xj^ when the 
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Figure 10: Allowed regions for the branching ratios of 7r+ — )■ e+ + (left), — )■ e+ + Ai"2,3 
(center), and — t- yU+ + A^2,3 (right) in terms of for the NH case. The regions within 
the red solid lines and the blue dotted lines are allowed for the case when r^vj a < 1 sec and 
TJV23 ''"bbn; respectively. The green horizontal lines represent the branching ratios when 
|0e/| = |©e/|uB (left and center) and |9^/| = |6^/|ub (right). 

conditions (ITSll and (flQil are satisfied, and hence the lower bounds on these branching ratios are 
also proportional to X'"^. For large values of the upper bounds on these branching ratios are 
smaller than those evaluated by using the experimental upper bound |Oe/|uB- This is because 
|9^/P must take its maximal value |6^/|ub ^'^ realize such a large X^, and then \Qei^^i\vB in 
addition to |6^/|ub forbids the case |6e/P = |0e/luB- other hand, the branching ratio 

of — 7- /i+ + A^2,3 cannot be so small since there is no cancellation in and also the larger 
value of 1 0^/1 is required for the shorter lifetime. Notice that |6^/p cannot be large as |B^/|ub 
for smaller values of X^. This is the reason why the upper bound on BR(K+ — > /i+ + X2_3) is 
smaller than that evaluated by |6^/p = |B^/|ub when X^ ^150. 

Therefore, we arrive at significant conclusions on the direct searches of X2.3 in the NH case. 
By the peak search of positrons in 7r+ or decays we would miss X2 and X3, since the model 
predicts too small branching ratios of these decays in some cases. Further, even if this is the 
case, it is possible to detect X2,3 in the yU"*" spectrum in the decays. Interestingly, the future 
experiments of BR(f^+ — + ^2,3) improved the sensitivity by two orders of magnitude can 
cover the whole parameter space when r^Vj 3 < 1 sec. 

In the IH shown in Fig. [TTl the regions of the branching ratios are different from 

those in the NH case. Especially, the branching ratios, vr e"*" + X2.3 and — )■ e"*" + ^2^3 
cannot be so small as in the NH case due to the absence of the cancellation in |0e/|. Further, 
it should be noted that the model predict the branching ratio of — )■ /i"*" + ^2,3 below the 
value with = IQ^^Iub the range of X^j. The reason for small values of X^^ is 

the same as the NH case. In addition, the ratio |6e/|/|0^j/| cannot be very small in the IH 
case, and then the stringent experimental bound |0e/|uB forbids being large as 1 0/^7-1 ub- 

This leads to the experimental signature of the IH case, i.e., the suppressed branching ratio of 
— )• /i"*" + ^2^3. In the IH case, therefore, the improvements of the experiments by one order 
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Figure 11: Allowed regions for the branching ratios of 7r+ — )■ e+ + (left), — e+ + Ai"2,3 
(center), and — )■ yU+ + N2^z (right) in terms of for the IH case. The regions with the red 
solid lines are allowed for the case when t^Tj 3 < 1 sec. The green horizontal lines represent the 
branching ratios when |9e/| = |6e/|uB (left and center) and = |6/x/|ub (right). 

of magnitude can cover the whole parameter range when t^^ 3 < 1 sec for all the three decay 
channels. 

Finally, we show in Figs. [12] and [13] the allowed ranges of the branching rations in terms of 
Mat for the NH and IH cases, respectively. Here we also impose the lifetime bound and vary 
the parameters of the model including X^. In the NH case, we find that the lower bounds on 
the branching ratios vr "'"—)■ e"*" + X2.3 and — > e"^ + ^2,3 are the same between two cases 
when r7V2 3 < ^ ^ec and rjvj.s < ''bbn- This is because such bounds are obtained by the largest 

as shown in Fig. [TU] It is also find that the upper bound for r^Va 3 < tbbn is slightly smaller 
than that for r^^ 3 < 1 sec. On the other hand, the upper bound on the branching ratio of 
— > /i+ + ^2,3 does not change by the lifetime bound. In contrast, the lower bound becomes 
much severer by the stringent lifetime bound tn^^ < tbbn- Therefore, the improvement in the 
measurement of BR(i^+ — )• /i+ + ^2,3) by a factor of 5 (by C(10^)) is sufficient to cover the 
allowed region with t^^ .^ < tbbn {'TN2.3 < 1 ^ec) for 60 MeV < < 130 MeV. In the IH case, 
for all three decay channels, the improvement in the measurements by one or two orders of 
magnitude allows to cover the allowed region for 60 MeV < M^v ^ 130 MeV. 

Before closing this section, we would like to remark the following two points. First, the 
properties of the mixing elements Qai obtained in Sec. [3] are independent of the choice of the 
masses M2,3. Therefore, the above arguments can be applied to sterile neutrinos heavier than 
pion mass. It should be stressed in particular that we have to pay a special attention to the 
processes of ^2,3 associated with electron/positron because of the cancellation in 9e/ for the 
NH case. Second, we would like to comment again that, although we consider the z/MSM in 
this paper, most of the results still hold for the general models of the seesaw mechanism with 
two right-handed neutrinos. 
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Figure 12: Allowed regions of the branching ratios of vr^ — )■ e"*" + iV2,3 (left), — )■ + N2,^ 
(center), and — i- /i"*" + (right) in terms of Mj^ for the NH case. The regions within 
the red solid lines and the blue dotted lines are allowed for the case when < 1 sec and 

Tn < T-BBN (USD. 
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Figure 13: Allowed regions of the branching ratios of vr^ — )■ e'^ + A'"2,3 (left), — )■ e"*" + A^2,3 
(center), and — >■ /x"*" + A''2,3 (right) in terms of Mjv for the IH case. The regions within the 
red solid lines are allowed with tat < 1 sec. 

6 Conclusions 

We have discussed mixing of active and sterile neutrinos in the i/MSM, paying special attention 
to the mixing elements, 0q,2 and 0q,3, of A'"2 and N^. In this model, these sterile neutrinos 
are responsible to generate the seesaw mass matrix of active neutrinos as well as the baryon 
asymmetry of the universe through the mechanism of neutrino oscillation. Since these mixing 
elements are crucial to determine the strength of interaction of A'"2,3, we have investigated the 
properties of Qa2 and 9^3 in detail. 

It has been shown that the parameter Imcj (or X^) is important to determine the overall 
scale of the Yukawa couplings of neutrinos. The couplings of N2,z scales as \Fai\ oc e'^™"^! and 
changes by orders of magnitude for |Imci;| ^ 1. We should stress again that the choice of Imo; 
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does not change the masses and mixing angles of active neutrinos. Thus, the mixing elements 
\Qai\ of A^2,3 can be larger as |Ima;| becomes large, being consistent with the oscillation data. 

We have then presented how the mixing elements depend on the neutrino parameters, 
namely masses, mixing angles and CP phases of active neutrinos. Interestingly, we have ob- 
served in the NH case that the leading terms of |6e2| and jOesI vanish at the same time and 
they becomes suppressed as Xj^ for the large X^^ region, if the conditions ( lT8l) and ( IT9l) are 
satisfied. This is our important result, which is illustrated in Fig. [TJ In order to realize this 
cancellation or the strong suppression in |9e2| and [Besl, the large value of ^13 is required. Thus, 
the experiments, e.g., Double Chooz, T2K, RENO, Daya Bay, and NOi^A in near future [65] 
will check whether it can happen in nature or not. The similar cancellation or suppression 
in |0e2| and {Qesl are potentially possible in the IH case. However, we have shown that the 
required conditions cannot be satisfied within the current neutrino data with 3a error. 

We have also discussed the Oz/2/3 decays in the i/MSM. Both contributions from active 
neutrinos and sterile neutrinos are fully taken into account. It has been found that quasi- 
degenerate sterile neutrinos N2 and generating the baryon asymmetry of the universe gives 
a significant effect when their masses are smaller than about 100 MeV which is the typical 
energy scale in the nuclear matrix elements of the decays. In fact, when their masses are well 
below 100 MeV, the effective neutrino mass in the Oz/2/3 decays vanishes in the z/MSM, as already 
pointed out in Ref. [50] by using the general seesaw model. Moreover, when their masses are 
comparable to lOOMeV, sterile neutrinos, N2 and A^s, give destructive contribution and then 
rries becomes much smaller than that from active neutrinos. In this case, therefore, the Oz/2/3 
decays are suppressed due to the existence of hght sterile neutrinos. Inversely speaking, the 
constraints from the Oi/2/3 decays on the parameters of the z/MSM become weaker in this mass 
region. On the other hand, when N2 and are much heavier than 100 MeV, the contributions 
to the Oz/2/3 decays from sterile neutrinos are negligible and rries is mainly comes from three 
active neutrinos as in the conventional seesaw model. 

Moreover, we have estimated the lifetimes of N2 and with M2.3 < m^r. It has been found 
that the lifetime can be shorter by taking large due to the enhancements of the mixing 
elements. As shown in Fig. [HI the present experimental limits on the mixing elements pose 
the lower bound on the lifetime as 10~^ sec at M^v — 120 MeV in the NH case and 10~^ sec 
at M^v — 120 MeV in the IH case. The obtained result is essential for the discussion of the 
cosmo logical constraints on the decays of A^2 and N3. Especially, we have pointed out that, even 
if one imposes the current cosmological bound from the BBN as tatj 3 < tbbn, there exists an 
allowed region for such light sterile neutrinos in the NH case. This is because of the cancellation 
in |0e/| obtained in the present analysis. We have then derived the predictions on neutrino 
parameters, namely, ^13 is large, ^ sm{6 + rj) c:^ —1, and ^ sini] < 0. When we take the weaker 
bound as 3 < 1 sec, there are the allowed regions for both NH and IH cases. 

Finally, we have investigated the direct search of A^2 and A3 with M2,3 < 171^^ in meson 
decays. Our study shows that, in the NH case, the experiments by using tt^ — )■ e"*" + A^2,3 
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and — > e+ + A^2,3 would miss A^2,3 when ^13 is large. This is because the cancellation or 
strong suppression in ©e/ can happen. In such a case, the peak search of the kaon decays 
+ A^2,3 is crucially important to find iV2,3 having sizable mixing elements Q^i. We 
should also comment that A^2,3 with Mat 55 MeV is hard to find in the NH case, if the 
cancellation in Q^i happens. This is because the peak search in — )■ /i"*" + A''2,3 is difficult 
for such a small mass region. On the other hand, in the IH case, all the three decay channels, 
7r+ -> e+H-A^2,3, e"'"+A^2,3, and A*^ + -^2,3, are promising for the future experiments 

finding A^2 and A^3. 
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A Mixing elements at 0{Xl) 

In this appendix, we shall present the expressions of the mixing elements |©a/P at the 0{X'^). 
In the NH case, we parameterize |©a/P at the 0{X'^) as 

|e.,n,5 = ^^.. (50) 
where A^, {a — e,fj,,T) are found as [r23 = (m2/m3)^/^] 

"--^'^ ^""-Q2^i3 + 2r23Csin(5 + r;)sin^i2tan^i3 + r^3 sin2^i2 , (51) 
^ /3 I „ t (.Qg g^^ (.Qg gjjj 2^23 sin 77 — sin 612 sin 2^13 sin^ ^23 sm{5 + rj) 

, (52) 



Ag = cos^ 6'i3 tan^ ^13 
Afj, = cos^ 6113 sin^ 623 + r23 C 



+^23 



Ar = COS^ 9i3 COS^ 6I23 - ^23 C 



+^23 



1 

cos^ 9i2 cos^ 6'23 + sin^ sin^ 613 sin^ ^^23 — sin 26'i2 sin 26*23 sin 613 cos 5 

COS Qx2 COS ^13 sin 2^23 sin 77 -|- sin Qyi sin 2^13 cos^ ^23 sin((5 -|- r\) 

^23 + 77 sin 2^12 sin 2^23 sin ^13 cos b . (53) 



cos^ Qx2 sin^ ^23 + sin^ Qyi sin^ 6*13 cos^ ^23 + \ 
On the other hand, in the IH case, we 



,21 _ ^2 
1^"-^! 1^2 ~ 4M7 



where 5^ [a = e, /x,r) are found as [ri2 = {m\jm2)^l'^\ 

Bg — cos^ 9i2 cos^ ^13 tan^ 9i2 — 2ri2^ sin 77 tan 612 + r^gj , 



(54) 



(55) 
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Bn = COS cos^ 6^23 + sin^ 6^12 sin Oh sin^ 6^23 — - sin 26*12 sin 26*23 sin ^13 cos 6 



[cos 26*12 cos 6 sin rj — sin 6 cos rj) sin 26*23 sin ^13 
+ sin 26*12 (cos^ 6*13 sin^ 6*23 + cos 26*23) sin r] 
cos^ 6*23 sin^ 6*12 + cos^ 6*12 sin^ 6*13 sin^ 623 + - sin 26*12 sin 26*23 sin 6*13 



(56) 



Br 



cos 6*12 sin^ 6*23 + sin^ 6*12 sin 6*13 cos^ 6*23 + - sin 26*12 sin 26*23 sin 6*13 cos 5 



( — cos 26*12 cos 6 sin 7] + sin 6 cos rj) sin 26*23 sin 6*13 
+ sin 26*12 (cos^ 6*13 cos^ 6*23 — cos 26*23) sinr] 



sin^ 6*23 sin^ 6*12 + cos^ 6*12 sin^ 6^13 cos^ 6*23 — - sin 26*12 sin 26*23 sin 6*13 



(57) 
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